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^—1 ' We consider the dynamics of a viscous cosmological fluid in the generalized Randall-Sundrum 

I model for an isotropic brane. To describe the dissipative effects we use the Israel-Hiscock-Stewart 

. full causal thermodynamic theory. In the limiting case of a stiff cosmological fluid with pressure 

04 ' equal to the energy density, the general solution of the field equations can be obtained in an exact 

f— I , parametric form for a cosmological fluid with constant bulk viscosity and with a bulk viscosity 

r-{ ■ coefficient proportional to the square root of the energy density, respectively. The obtained solutions 

' describe generally non-inflationary brane worlds, starting from a singular state. During this phase 

QQ , of evolution the comoving entropy of the Universe is an increasing function of time, and thus a large 

' amount of entropy is created in the brane world due to viscous dissipative processes. 

T— I ! 

> ■ I. INTRODUCTION 

m ; 

\ The idea [|l|,H that our four-dimensional Universe might be a three-brane embedded in a higher dimensional space- 
\^ • time has recently attracted much attention. According to the brane-world scenario, the physical fields in our four- 
<^ ' dimensional space-time, which are assumed to arise as fluctuations of branes in string theories, are confined to the three 
: brane, while gravity can freely propagate in the bulk space-time, with the gravitational self-couplings not significantly 
modified. This model arisen from the study of a single 3-brane embedded in five dimensions, with the 5D metric 
[ given by ds"^ = e~^^y^ri^i,dx^dx^ + dy^, which can produce a large hierarchy between the scale of particle physics 
■ and gravity due to the appearance of the warp factor. Even if the fifth dimension is uncompactified, standard AD 
O il gravity is reproduced on the brane. Hence this model allows the presence of large or even infinite non-compact extra 
^ ' dimensions and our brane is identified to a domain wall in a 5-dimensional anti-de Sitter space-time. 

The Randall-Sundrum model was inspired by superstring theory. The ten-dimensional Eg x Eg heterotic string 
theory, which contains the standard model of elementary particle, could be a promising candidate for the description 
' of the real Universe. This theory is connected with an eleven-dimensional theory, M-theory, compactified on the 
rN : orbifold X 5*7^2 [|. In this model we have two separated ten-dimensional manifolds. For a recent review of 
• dynamics and geometry of brane Universes see [|j . 

The static Randall-Sundrum solution has been extended to time-dependent solutions and their cosmological prop- 
erties have been extensively studied ^ - jlj]. In one of the first cosmological applications of this scenario it was 
pointed out that a model with a non-compact fifth dimension is potentially viable, while the scenario which might 
solve the hierarchy problem predicts a contracting Universe, leading to a variety of cosmological problems [0 . By 
adding cosmological constants to the brane and bulk, the problem of the correct behavior of the Hubble parameter 
on the brane has been solved by Cline, Grojean and Servant ]T8| . As a result one also obtains normal expansion 
during nucleosynthesis, but faster than normal expansion in the very early Universe. The creation of a spherically 
symmetric brane-world in AdS bulk has been considered, from a quantum cosmological point of view, with the use of 
the Wheeler-de Witt equation, by Anchordoqui, Nunez and Olsen flst] . 

The effective gravitational field equations on the brane world, in which all the matter forces except gravity are 
confined on the 3-brane in a 5-dimensional space-time with Z2-symmetry have been obtained, by using an elegant 
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geometric approach, by Shiromizu, Maeda and Sasaki |20 2^. The correct signature for gravity is provided by the 
brane with positive tension. If the bulk space-time is exactly anti-de Sitter, generically the matter on the brane is 
required to be spatially homogeneous. The electric part of the 5-dimensional Weyl tensor Ejj gives the leading order 
corrections to the conventional Einstein equations on the brane. The four-dimensional field equations for the induced 
metric and scalar field on the world-volume of a 3-brane in the five-dimensional bulk with Einstein gravity plus a 
self- interacting scalar field have been derived by Maeda and Wands ||2^. The effective four-dimcnsional Einstein 
equations include terms due to scalar fields and gravitational waves in the bulk. 

The linearized perturbation equations in the generalized Randall-Sundrum model have been obtained, by using 
the covariant nonlinear dynamical equations for the gravitational and matter fields on the brane, by Maartens [ p3[ . 
The nonlocal energy density determines the tidal acceleration in the off-brane direction and can oppose singularity 
formation via the generalized Raychaudhuri equation. Isotropy of the cosmic microwave background may no longer 
guarantee a Friedmann-Robcrtson- Walker geometry. Vorticity on the brane decays as in general relativity, but nonlocal 
bulk effects can source the gravitomagnetic field, so that vector perturbations can also be generated in the absence of 
vorticity. 

A systematic analysis, using dynamical systems techniques, of the qualitative behavior of the Friedmann-Robertson- 
Walker (FRW), Bianchi type I and V cosmological models in the Randall-Sundrum brane world scenario, with matter 
on the brane obeying a barotropic equation of state p = — l)p, has been realized by Campos and Sopuerta p^ , ^ . 
In particular, they constructed the state spaces for these models and discussed what new critical points appear, the 
occurrence of bifurcations and the dynamics of the anisotropy for both a vanishing and non- vanishing Weyl tensor in 
the bulk. 

All these investigations of brane cosmological models have been performed under the simplifying assumption of a 
perfect cosmological fluid. But in many cosmological situations an idealized fluid model of matter is inappropriate, 
especially in the case of matter at very high densities and pressures. Such possible are the relativistic transport of 
photons, mixtures of cosmic elementary particles, evolution of cosmic strings due to their interaction with each other 
and surrounding matter, classical description of the (quantum) particle production phase, interaction between matter 
and radiation, quark and gluon plasma viscosity etc. From a physical point of view the inclusion of dissipative terms 
in the energy-momentum tensor of the cosmological fluid seems to be the best motivated generalization of the matter 
term of the gravitational field equations. 

The first attempts at creating a theory of relativistic fluids were those of Eckart j2^ and Landau and Lifshitz 
|p7| . These theories are now known to be pathological in several respects. Regardless of the choice of equation of 
state, all equilibrium states in these theories are unstable and in addition signals may be propagated through the 
fluid at velocities exceeding the speed of light. These problems arise due to the first order nature of the theory, that 
is, it considers only first-order deviations from the equilibrium leading to parabolic differential equations, hence to 
infinite speeds of propagation for heat flow and viscosity, in contradiction with the principle of causality. Conventional 
theory is thus applicable only to phenomena which are quasi-stationary, i.e. slowly varying on space and time scales 
characterized by mean free path and mean collision time. 

A relativistic second-order theory was found by Israel and developed by Israel and Stewart into what is 
called "transient" or "extended" irreversible thermodynamics. In this model deviations from equilibrium (bulk stress, 
heat flow and shear stress) are treated as independent dynamical variables, leading to a total of 14 dynamical fluid 
variables to be determined. However, Hiscock and Lindblom |30j and Hiscock and Salmonson have shown that 
most versions of the causal second order theories omit certain divergence terms. The truncated causal thermodynamics 
of bulk viscosity leads to pathological behavior in the late Universe while the solutions of the full causal theory are Q : 
a) for stable fluid configurations the dissipative signals propagate causally, b) unlike in Eckart-type's theories, there 
is no generic short wave-length secular instability and c) even for rotating fluids, the perturbations have a well-posed 
initial value problem. For general reviews on causal thermodynamics and its role in relativity see and p^ . 

Causal bulk viscous thermodynamics has been extensively used for describing the dynamics and evolution of the 
early Universe or in an astrophysical context. But due to the complicated character of the evolution equations, very 
few exact cosmological solutions of the gravitational field equations are known in the framework of the full causal 
theory. For a homogeneous Universe filled with a full causal viscous fluid source obeying the relation ^ ^ P^^^, 
exact general solutions of the field equations have been obtained in |^5|-^. In this case the evolution of the bulk 
viscous cosmological model can be reduced to a Painleve-Ince type differential equation, whose invariant form can be 
reduced, by means of nonlocal transformations, to a linear inhomogeneous ordinary second-order differential equation 
with constant coefficients It has also been proposed that causal bulk viscous thermodynamics can model on a 
phenomenological level matter creation in the early Universe |^,^. Exact causal viscous cosmologies with ^ ~ 
have been obtained in |Q . 

Because of technical reasons, most investigations of dissipative causal cosmologies have assumed FRW symmetry 
(i.e. homogeneity and isotropy) or small perturbations around it pl[ |. The Einstein field equations for homogeneous 
models with dissipative fluids can be decoupled and therefore are reduced to an autonomous system of first order 
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ordinary differential equations, which can be analyzed qualitatively | |42| , [43| . 

It is the purpose of the present paper to investigate the effects of the bulk viscosity of the cosmological matter 
fluid on the dynamics of the brane world. Since the effects of the extra-dimensions and also the viscous effects are 
more important at high matter densities, we restrict our analysis to the extreme case of a stiff (Zeldovich type) 
cosmological fluid, with pressure equal to the energy density. Hence the most important contribution to the energy 
density of the matter comes from the quadratic term in density, and during this period the energy density of matter 
is proportional to the Hubble parameter, in opposition to the standard general relativistic case with energy density 
proportional to the square of the Hubble parameter. In this case, and by assuming that the bulk viscosity coefficient 
and the temperature dependence of the cosmic fluid on the energy density are given by a simple power laws, the fleld 
equations can be solved exactly for several explicit functional forms of the viscosity coefficient. 

The present paper is organized as follows. The fleld equations on the brane describing the evolution of a viscous 
cosmological fluid arc written down in Section II. In Section HI we present the general solution of the field equations 
for a constant bulk viscosity and a bulk viscosity coefficient proportional to the square root of the energy density. 
The study of the stability of the equilibrium points of the dynamical system associated to the evolution of the viscous 
cosmological fluid is performed in Section IV. In Section V we discuss and conclude our results. 



II. DISSIPATIVE COSMOLOGICAL FLUIDS ON THE BRANE 



In the hD space-time the brane-world is located as Y{X^) — 0, where , I = 0,1,2,3,4 are 5-dimensional 
coordinates. The effective action in flve dimensions is 

-'^ D All/" / - ( 1 j,^± \ I r matter A 



j d^XV^, (^^i?5 - As j + J^_^ ct'x^g ^-j^K^ - A + L™J , (1) 

with fc| = SttGs the 5-dimensional gravitational coupling constant and where x^, ^ = 0,1,2,3 are the induced 4- 
dimensional brane world coordinates. i?5 is the 5D intrinsic curvature in the bulk and is the intrinsic curvature 
on either side of the brane. 

On the 5-dimensional space-time (the bulk), with the negative vacuum energy A5 and brane energy-momentum as 
source of the gravitational field, the Einstein field equations are given by 

Gij = klTij, Tij = -A.gij + S{Y) [-Xgu + T.^""] , (2) 

In this space-time a brane is a fixed point of the Z2 symmetry. In the following capital Latin indices run in the range 
0, 4 while Greek indices take the values 0, 3. 

Assuming a metric of the form ds^ — {ninj + gjj)dx^dx'^, with njdx^ = dx the unit normal to the x = const, 
hypersurfaces and gu the induced metric on x = const, hypersurfaces, the effective four-dimensional gravitational 
equations on the brane (which are the consequence of the Gauss-Codazzi equations) take the form ]20| , ^ : 

Gf_ii, — —Ag^i, + k^T^i, -|- k'^S^u — E^^^ (3) 

where S^j^i, is the local quadratic energy-momentum correction 

V = ^TT,, - -T.'^T^c + (3T"'5t„^ - T^) , (4) 

and E^^ is the nonlocal effect from the bulk free gravitational filed, transmitted projection of the bulk Weyl tensor 

ClAJB 

Eij ^CiAjBn^n'', Ejj^E^J'^S'^j as x ^ 0. (5) 
The four-dimensional cosmological constant. A, and the coupling constant, ki, are given by 



A = f^A. + ^J, ki^^. (6) 

The Einstein equation in the bulk, Codazzi equation, also implies the conservation of the energy momentum tensor 
of the matter on the brane, 

D^T^" = 0. (7) 
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Moreover, the contracted Bianchi identities on the brane imply that the projected Weyl tensor should obey the 
constraint 



(8) 



Finally, the equations (|^, ^ and give the complete set field equations for the brane gravitational field. 

For any matter fields (scalar field, perfect fiuids, kinetic gases, dissipative fiuids etc.) the general form of the brane 
energy-momentum tensor can be covariantly given as 



7)11/ = pu^u^ + ph 



(9) 



The decomposition is irreducible for any chosen 4- velocity it^. Here p and p are the energy density and isotropic 
pressure, and h^^ = + u^u^ projects orthogonal to u^. The energy flux obeys = <z<^>, and the anisotropic 



stress obeys tt 



'<fj.u>7 where angular brackets denote the projected, symmetric and tracefree part: 



W, 



,/3 
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(10) 



The symmetric properties of E^i, imply that in general we can decompose it irreducibly with respect to a chosen 
4-velocity field m'' as 



E., 



U 



(11) 



where k = k^/k^. In Eq. ( |ll| ) U is the effective nonlocal energy density of the brane arising from the free gravitational 
field in the bulk, Vfj,i, is the nonlocal anisotropic stress carying Coulomb, gravito-magnetic and gravitational wave 
effects from the bulk, while Q is the effective nonlocal energy flux on the brane. 

The effect of the bulk viscosity of the cosmological fluid can be considered by adding to the usual thermodynamic 
pressure p the bulk viscous pressure 11 and formally substituting the pressure terms in the energy-momentum tensor 
by PeS = p + n. The particle flow vector N'^ is given by N'^ — nu'^\ where n > is the particle number density. 

In the framework of causal thermodynamics and limiting ourselves to second-order deviations from equilibrium, the 
entropy flow vector takes the form 



Tip 



(12) 



where s is the entropy per particle, t the relaxation time, T the temperature and ^ is the bulk viscosity coefficient. 

We consider that the heat transfer is zero, 9^ = in (g), and also a vanishing effective nonlocal anisotropic stress 
and energy fiux, 7^^,^ = = Q,,. Then the matter corrections are given by 
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The line element of a fiat Robertson- Walker metric is given by 

ds^ = -dt'^ + a^{t) {dx^ + dy^ + dz^) . 



(13) 



(14) 



We also assume that the thermodynamic pressure p of the cosmological fluid obeys a linear barotropic equation of 
state p — {'y — l)p, 7 = const, and 1 < 7 < 2. Under these assumptions, the field equations and the conservation 
equations for the Bianchi type I brane gravitational field take the form 



= A 



9 /b? 9 6Z^ 



2H 



3H' = A- kl [(7 - l)p + n] - ^ [(27 - l)p2 



k^X 



p + 3jHp = -3i?n, 
U + AHU = 0, 

where the Hubble parameter H is defined as H — d/a. N'^-^ 
h + 3Hn 0. 



(15) 

(16) 

(17) 
(18) 



leads to the particle number conservation equation 
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The causal evolution equation for the bulk viscous pressure 11 is given by 

rtl + Il = -3^H-^Tn(^3H+^-^-^y (19) 

Eq. ( p^ arises as the simplest way (linear in H) to satisfy the iJ-theorem (i.e., for the entropy production to be non- 
negative, 5''^;;^ = lP/{£,T) > 0). The Israel-Stewart theory is derived under the assumption that the thermodynamic 
state of the fluid is close to equilibrium, which means that the non-equilibrium bulk viscous pressure should be small 
when compared to the local equilibrium pressure, that is |n| < p. 

The growth of the total comoving entropy S(t) over a proper time interval {to, t) is given by [p3 



S(i)-S(to)=-T^ f^^dt. (20) 



An important observational quantity is the deceleration parameter q — dH~^ /dt — 1. The sign of the deceleration 
parameter indicates whether the model inflates or not. The positive sign of q corresponds to "standard" decelerating 
models whereas the negative sign indicates inflation. 

Since the effects of the extra-dimensions are important at very high densities, when the cosmological fluid behaves 
like a Zeldovich fluid with p = p (7 = 2), as are also the dissipative effects, we consider only the physical situation in 
which the quadratic term dominates in the energy equation (p^). Therefore during the early period of evolution the 
energy density of the Universe is given approximately by p w 3po-ff , with 3po = a/ 6A/fc|. 

In order to close the system of equations (p^-(p^ we have to specify T, r and ^. 

First, following |44j we suppose that the relation t = ^/p holds in order to guarantee that the propagation velocity 
of bulk viscous perturbations, i.e. the nonadiabatic contribution to the speed of sound in a dissipative fluid without 
heat flux or shear viscosity does not exceed the speed of light. An analysis of the relativistic kinetic equation for 
some simple cases given by Belinskii and Khalatnikov |45j, Belinskii, Nikomarov and Khalatnikov Q] and Murphy 
p6[ has shown that in the asymptotic regions of small and large values of the energy density, the viscosity coefficients 
can be approximated by power functions of the energy density with definite requirements on the exponents of these 
functions. For small values of the energy density it is reasonable to consider large exponents, equal in the extreme 
case to one. For large p the power of the bulk viscosity coefficient should be considered smaller (or equal) to 1/2. 

Therefore we assume the following simple phenomenological laws for the bulk viscosity coefficient, temperature and 
relaxation time: 

^^ap'^CoH^ T^Pp^^ToH\ r=^ = ^^, (21) 

P To 

where s>0, r>0, q;>0 and /3 > are constants and ^0 — Q!(3po)* and tq = £,o^- 

In the context of irreversible thermodynamics, p, p, T and the number density n are equilibrium magnitudes which 
are generally related by equations of state of the form p = p(T, n) and p — p{T,n). From the requirement that the 
entropy is a state function we obtain the equation 

dp\ p + P T f dp 



» I \ ' ' (22) 

on / rp n n \oi 

which imposes the constraint r = (7 — l)/7- Hence for a Zeldovich fluid we have r = 1/2. 

With these assumptions the bulk viscous pressure 11 can be obtained from Eq. ( p7[ ) in the form 

n = -po l^l + eiJ^ (23) 

The bulk viscous evolution equation ( |l9[ ) can be written as 

U+-Il = -3pH-ln(3H-l^) , (24) 



r 2 V 2p 

and for a stiff cosmological fluid on the brane takes the form 

ii 7 



^---^ + {3 + toH-'^)H + QtoH^-' = Q. (25) 
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By means of the substitution H = y Eq. ( p5[ ) takes the form 

j/ + (3 + roy'^/') y-'lH - ^roy(4^-^)/3 = 0. (26) 

By taking u — y and denoting v — Eq. (|2^) can be transformed to a second type Abel nonhnear first order 
differential equation: 



dv 
dy 

By introducing a new variable 77 — Sj/^^^'^, Eq. (E7h becomes 



^ - (3 + T,y'^/^) y-''\' + ^roy(^^-^)/3z;3 = 0. (27) 



^ + (3 + S^Vory-^^) - ^ro^^-^^^^^ = 0. (28) 



III. BRANE EVOLUTION OF DISSIPATIVE STIFF COSMOLOGICAL FLUIDS 

In the previous Section we have formulated the basic equations describing the dynamics of a dissipative stiff 
cosmological fluid on the brane. We have considered only the extreme case of very high densities, when the main 
contribution to the energy of the matter is given by the quadratic term in the energy-momentum tensor, due to the 
form of the Gauss-Codazzi equations, and which leads to major changes in the dynamics of the Universe. In this case 
the basic equation describing the evolution of the Universe can be reduced to an Abel type equation (|2^). 

It is the purpose of the present Section to consider some exact classes of solutions of Eq. ( |28| ) , corresponding to 
some particular values of the constant s. 

As a first case we assume that the bulk viscosity coefficient ^ is a constant, ^ = ^0 = const., corresponding to the 
choice s = in the equation of state of the bulk viscosity coefhcient. For s = the temperature and the relaxation 
time are functions of density, according to the equations of state (|2l|). 

For s = the evolution equation (Eq) of the bulk viscous pressure takes the form 

^ + (3 + ro)z;2-^w3^0. (29) 
drj 2 

By introducing a new variable 77' = (3 + ro)?y and denoting b = 3to/2(3 + tq)^, Eq. ( p9| ) takes the form 

^+v^- bTj'v^ = 0. (30) 
drj' 



By taking v — w/rj' , Eq. (30) is transformed into 

dw w 



d,' ,7(^-^ -- + !)■ (31) 



Hence the general solution of Eq. (|3y) is given by 



^/-g / =e/W, (32) 
V bw^ — w + 1 



where C > is an arbitrary constant of integration. 



and 



f{w) = In ^ , if 6 < -, 34 

^ ^ 2^/A \2bw-l + ^J '4' ^ ' 

/H = -:t%' if b=\, (35) 



w ■ 



'^^^^ " arctan ^^^^=i j , if & > ^, (36) 
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where we denoted A = 1 — 46. 

Therefore the general solution of the field equations can be expressed in the following exact parametric form, with 
9 = l/w taken as parameter: 



tie) - to - Co I [e^ -9 + b)e-'^f'-^Ue, 



aie)=aoexp[-2CoHofie)] 
q{e) = Ae/CoHo - 1, 



11(9) = -2po- 



»4/(e) 



:{3Ho-2e/Co)., 



Kb To 



?2 - 6* + 5)2 

^ Galpo^ f (SCoHo - 26) exp[2/(0)(l - 3Coi/o) 



(612 _ 6) + 6)2 



■d9. 



(37) 

(38) 

(39) 
(40) 

(41) 

(42) 
(43) 



where Hq — [C/3(3 + to)]"', Co — 3'*(3 + to)'^/C^ and to, qq and Uq are constants of integration. In the new variable 9 
the function / is given by f{0) = -2'^ J{9^ -9 + b)-^d9. 

The thermodynamic consistency of the model can be studied from the ratio of the bulk viscous and thermodynamic 
pressure, which is given by 



|5-g| 



(44) 



The second case we analyze corresponds to the extreme limit of very high densities when s — 1/2. Then Eq. (|2S 
takes the form 



dv 
drj 



+ (3 + QTorf 



27 



Introducing two new functions A{ri) = -~2ri'^/9TQ + 2/3 and B{r]) = — 2?7/27to allows to rewrite Eq. 
general form 



dv 



dr] B{ri) 



d A{rj) 



drj B{rj) 



By introducing a new variable 



1 A{r^) 
a = 



Eq. (Bq) can be written in the general form 



or, equivalently, 



dr] 
da 



V B{r]) 



^ ^ B{r,)a + Air,), 
da 



2^2 2 
"9n)^ ~ 27^^^*^^ 3' 



(45) 
in the 

(46) 

(47) 

(48) 
(49) 



Hence we have transformed the initial Abel type equation into a Riccati equation. A particular solution of Eq. (|4£ 
is given by 



rj = 9rocrA((7), 



(50) 



and therefore the general solution of Eq. ( [49|) is 
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rj{a) = 9rocrA(cr) 



A2(a)e-"'/27-o 



Ci + g|^/A2(a)e-V27rod^' 



(51) 



where A(cr) = (cr^ + 27to/2)-\ 

Hence we obtain the general solution of the gravitational field equations on the brane for a Zeldovich causal bulk 
viscous fluid, with bulk viscosity coefficient proportional to the square root of the density, in the following exact 
parametric form, with a taken as parameter: 



t{a) ^ to = 6t-^ J ri-^{a)da, 



a{(7) — flo exp 



Act 



Ma) 
27ro 
36ro , 



77(0-) 7?2((t) 



da 
11. 



81 

= E(ao) - 
U (a) — Uq exp 



8lToroA;B 
8 

"27^ 



r]ia)[3r]^{a) + 2arj{a) - IStq] exp 



7]{a)da 



9to 



da 



da. 



(52) 
(53) 
(54) 
(55) 
(56) 
(57) 
(58) 



IV. STABILITY ANALYSIS OF THE EQUILIBRIUM POINTS OF THE VISCOUS COSMOLOGICAL 

FLUID 



The general evolution equation of the bulk viscous cosmological fluid on the brane is given by Eq. (^5|). From 
mathematical point of view it is a second order non-linear differential equation of the form H + R{H, H) = 0, with 
R{H, H) = -IH'^/AH + {iH + tqH^-^)H + &tqH^-\ Therefore i?(0, 0) = Yann^x^o R{H, X) = and R{H, 0) 7^ 
for H ^Q. 

In order to study the stability of the equilibrium points of the evolution equation of the viscous cosmological fluid 
on the brane, Eq. (p5|), we shall rewrite it in the form of an autonomous dynamical system, by introducing a new 
variable X = H: 

dH , , 

dX 7 

^ = L--{3H + roH'-nX~6ToH'-^. (60) 

The critical points of this dynamical system are given hy H = X = 0. They corresponds to a Minkowskian space- 
time (a = const. = 1) and to a de Sitter inflationary phase, with a — exp{HQt), Hq = const.. The system has no other 
critical points besides the origin. 

The Lyapunov function V{H,X) associated to this system can be chosen Q as V{H, X) = X^/2 + j" R{s, 0)ds 
and is given by 

F(iJ,X) = ix2 + ^ij4-. (61) 



The Lyapunov function (|l]) has the properties V{0,0) = and dV/dt = 7X^/AH - {3H + toH^-')X^. According 
to the standard theory of this type of differential equations the equilibrium state {H = 0, X = 0) is stable if the 
conditions 

HR{H, 0) = 6toH'^-' > 0, for H^O, (62) 

and 
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dR{H,X) 
dX 



7X 
2H 



+ 1) + 3 + tqH- 



> 0, 



(63) 



Moreover, if the condition dR{H, X)/dX > is satisfied for 
M. The equilibrium state is unstable if dR{H,X)/dX < 



holds, where the deceleration parameter q = —X/H^ — 1 
HX 0, the equilibrium state is asymptotically stable 
for HX ^0. ~ 

The stability criteria of the critical point can be formulated in terms of some conditions imposed on the deceleration 
parameter. In the limit of large time, H ^ and the term > dominates in the expression of dR{H,X)/dX, 
making it obviously non-negative. In the small time limit, iJ — > oo and the condition of the stability of the critical 
point is 7{q + l)/2 + 3 > 0, or q > —13/7. If g > —13/7 the critical point is also asymptotically stable. On the other 
hand for q < —13/7 the equilibrium point is unstable. 



V. DISCUSSIONS AND FINAL REMARKS 



In the present paper we have considered the evolution of a causal viscous dissipative cosmological fluid in the brane 
world scenario. As only source of dissipation we have considered the bulk viscosity of the matter on the brane. The 
most important differences to standard general relativity are expected to occur in the limit of extremely high densities, 
when the fluid obeys a Zeldovich (stiff) equation of state p = p. In this case the Friedmann equation is modified 
due to the presence of the terms from extra dimensions, quadratic in the energy density, which dominates the other 
terms in the energy momentum tensor, leading to an energy density of the cosmic fluid proportional to the Hubble 
parameter. 

By assuming the usual equations of state for bulk viscosity, temperature and relaxation time, the field equations can 
be solved exactly for some specific choices of the constant s describing the bulk viscosity coefRcient- energy density 
functional relation. 

For s = 0, case corresponding to a constant bulk viscosity coefficient ^ = const., the general solution of the field 
equations for the viscous fluid on the brane world is given by Eqs. (p7|)-(^3|). Since the bulk viscous pressure 11 must 
be negative, 11 < 0, it follows that in order to satisfy this condition the range of the parameter must be restricted 
to 61 e [0,2/3]. 

In the limit of small times, we have 6 —^ and one obtains the following equations describing the evolution of the 
viscous cosmological fluid on the brane: 

jj^^-2CoHot^ p ^ p t^'^CoHot ^ 9 = -1, 

n Z^-i-2CoHo^ T,{t) ^ll{to) +t^^^-^°"°^/'^. (64) 

The variations of the scale factor, Hubble parameter, deceleration parameter, bulk viscous pressure and comoving 
entropy of the constant bulk viscosity coefficient dissipative cosmological fluid confined on the brane are presented in 
Figs. 1-5. 
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FIG. 1. Variation as a function of time of the Hubble parameter H of the brane Universe with confined dissipative cosmological 
fluid with constant coefficient of bulk viscosity for different values of the parameter 6: 6 = 1/6 (sofid curve), 6 = 1/4 (dotted 
curve) and 6=1/2 (dashed curve). The values of the constants Ho, Co and to are different for each curve. 
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FIG. 2. Time evolution of the scale factor a of the branc Universe with confined dissipative cosmological fluid with constant 
coofficicut of bulk viscosity for different values of the parameter b: & = 1/6 (solid curve), 6 = 1/4 (dotted curve) and 6 = 1/2 
(dashed curve). The values of the constants ao, Co and to are different for each curve. 




FIG. 3. Dynamics of the deceleration parameter q of the brane Universe with confined dissipative cosmological fluid with 
constant coefficient of bulk viscosity for different values of the parameter 6: 6 = 1/6 (solid curve), 6=1/4 (dotted curve) 
and 6 = 1/2 (dashed curve) . The value of the constant to is different for each curve and the constants Co and Ho have been 
normalized so that HoCo = 1. 
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FIG. 4. Variation as a function of time of the bulk viscous pressure 11 of the brane Universe with conflned dissipative 

cosmological fluid with constant coefficient of bulk viscosity for different values of the parameter b: b = 1/6 (solid curve), 
6=1/4 (dotted curve) and 6=1/2 (dashed curve). The values of the constants Ho, Co and to are different for each curve. 
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FIG. 5. Time evolution of the comoving entropy E of the brane Universe with confined dissipative cosmological fiuid with 
constant coefficient of bulk viscosity for different values of the parameter b: b = 1/6 (solid curve), b — 1/4 (dotted curve) and 
b = 1/2 (dashed curve). The values of the constants Ho, Co, to and ao are different for each curve. 

The evolution of the viscous brane Universe starts from a singular state, with infinite values of the energy density 
and pressure and zero scale factor. The initial evolution is inflationary, with a negative deceleration parameter. But 
in the large time limit the dynamics becomes non-inflationary, with the deceleration parameter q > 0. Therefore 
the inclusion of viscous effects during the period when the quadratic term in energy density (coming from extra- 
dimensions) dominates the dynamics of the space time, provides an effective mechanism for the "graceful exit" of the 
brane world from the initial inflationary phase to a non-inflationary era. Due to the dissipative effects the entropy on 
the brane is increasing in time and a large amount of entropy is produced in the early stages of the evolution of the 
brane Universe. The nonlocal energy density on the brane, U, is a decreasing function of time, so the effects of the 
gravitational field on the bulk become rapidly negligeable. But the criterion of the thermodynamic consistency of the 
model, I = |n/p| < 1, is not generally satisfied in this model, as can be easily seen from Eq. (^). All the inflationary 
states clearly contradicts the condition. On the other hand it is possible to find some particular of the parameter b 
leading to thermodynamic consistency during the noninflationary phase. 

The time variations of the Hubble parameter, scale factor, deceleration parameter, bulk viscous pressure and 
comoving entropy for the brane Universe with a dissipative cosmological fluid with the bulk viscosity coefficient 
proportional to the square root of the energy density (s = 1/2) are represented in Figs. 6-10. 
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FIG. 6. Variation as a function of time of the Hubble parameter H of the brane Universe with confined dissipative cosmological 
fluid with the coefficient of bulk viscosity proportional to the square root of the energy density (s = 1/2) for different numerical 
values of the integration constant Ci: Ci = 0.58 (sofid curve), Ci = 0.62 (dotted curve) and Ci = 0.66 (dashed curve) {to = 0). 
The constant to has been normalized so that tq — 1. 
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FIG. 7. Time evolution of the scale factor a of the brane Universe with confined dissipative cosmological fluid with the 
coefficient of bulk viscosity proportional to the square root of the energy density (s — 1/2) for different numerical values of 
the integration constant Ci: Ci = 0.58 (solid curve), Ci = 0.62 (dotted curve) and Ci = 0.66 (dashed curve) (to = 0). The 
constant To has been normalized so that To = 1. 
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FIG. 8. Dynamics of the deceleration parameter q of the brane Universe with confined dissipative cosmological fiuid with 
the coefficient of bulk viscosity proportional to the square root of the energy density (s = 1/2) for different numerical values 
of the integration constant Ci: Ci = 0.58 (solid curve), Ci = 0.62 (dotted curve) and Ci = 0.66 (dashed curve) {to = 0). The 
constant To has been normalized so that ro = 1. The deceleration parameter satisfies the condition g > for all times. 
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FIG. 9. Variation as a function of time of the bulk viscous pressure 11 of the brane Universe with confined dissipative 

cosmological fluid with the coefficient of bulk viscosity proportional to the square root of the energy density (s = 1/2) for 
different numerical values of the integration constant Ci: Ci = 0.58 (solid curve), Ci = 0.62 (dotted curve) and Ci = 0.66 
(dashed curve) (to = 0). The constant To has been normalized so that ro = 1. As required by the model, the bulk viscous 
pressure is negative for all times. 
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FIG. 10. Dynamics of the comoving entropy E of the brane Universe with confined dissipative cosmological fluid with the 
coefficient of bulk viscosity proportional to the square root of the energy density (s — 1/2) for different numerical values of 
the integration constant Ci: Ci — 0.58 (solid curve), Ci = 0.62 (dotted curve) and Ci = 0.66 (dashed curve) (to = 0). The 
constant to has been normalized so that tq — 1. 



The behavior of the Universe depends on the numerical values of the arbitrary integration constant Ci and of the 
constant tq. For the chosen numerical values of these parameters the Universe generally starts from a singular state, 
with zero and infinite values of the scale factor and energy density, respectively. In order the model represents a 
dissipative fluid, with negative bulk viscous pressure, the parameter a must satisfy the condition 2a < ISro/rj^a) — 
3rj{(7). In opposition to the constant bulk viscosity case, the evolution is non-inflationary for all times. Due to viscous 
dissipative effects a large amount of comoving entropy is created on the brane and the entropy of the Universe is 
increasing due to viscous dissipation. 

Fig. 11 presents the time variation of the ratio I of the bulk viscous and thermodynamic pressures, respectively. 
For all times |n/p| < 1 and hence in this model the thermodynamical consistency condition of the smallness of the 
bulk viscous pressure is satisfied for all times and for a large class of admissible values of the integration constant Ci 
and of Tq. 
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FIG. 11. Time variation of the absolute value / of the ratio of the bulk viscous and thermodynamic pressures, I — |n/p| of 
the brane Universe with confined dissipative cosmological fluid with the coefficient of bulk viscosity proportional to the square 
root of the energy density (s — 1/2) for different numerical values of the integration constant Ci: Ci — 0.58 (solid curve), 
Ci = 0.62 (dotted curve) and Ci = 0.66 (dashed curve) (to = 0). The constant tq has been normalized so that tq — 1. For all 
times for the chosen set of parameters the ratio of the pressures is smaller than 1. 



The effective nonlocal energy on the brane, U, is tending rapidly to zero in the large time limit. 

The general exact solution of the gravitational field equations for a homogeneous flat FRW Universe filled with 
a causal bulk viscous fluid with the bulk viscosity coefficient proportional to the Hubble function, f ~ p^/^ ^ H 
has been obtained, in the framework of the standard general relativity (GR), in |^^. The solution of the field 
equations can in this case also be represented in an exact parametric form. There are major differences between 
the general evolution of the dissipative cosmological fiuids in the brane and GR models. In GR the dynamics of the 
cosmological fluid described by the general solution is inflationary for all times, with the thermodynamic consistency 
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condition \Tl/p\ < 1 violated during the entire expansionary evolution period. But in the brane model the evolution 
is non-inflationary, with the consistency condition satisfied, at least for a specific range of values of the parameters 
To and Ci, which are unknown for a realistic physical situation (of course for some particular numerical values of 
these parameters inflationary initial behavior or increasing energy density can also be obtained). During the general 
relativistic inflationary period the density, temperature, bulk viscosity coefficient and comoving entropy are rapidly 
increasing functions of time. In fact the general solution of the GR field equations describe a transition between two 
Minkowskian space-times connected by an inflationary period. For some particular values of the parameters one can 
also obtain general relativistic non-inflationary solutions \ ^7\ . 

The consideration of viscous dissipative efects in the brane and general relativistic models in the extreme limit 
of very high densities could be a useful way to differentiate between the two cosmological scenarios. The different 
behavior of the energy density of cosmic matter (p ~ iJ in brane and p ~ in GR, respectively), leads, via the 
bulk viscous pressure evolution equation, to differences in the dynamics of the very early Universe, which perhaps can 
serve as a tool for testing the viability of brane model cosmology. 
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